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' Abstract 

We show that the Adomian decomposition method, the time-series 
expansion, the homotopy-perturbation method, and the variational- 
iteration method completely fail to provide a reasonable description 
■ of the dynamics of the simplest prey-predator system. 

OO ' 

In is well known that a reasonable description of the dynamics of a 
nonlinear system by means of a perturbation method is a difficult enter- 
prize [1]. Recently, several approaches have been proposed for the treatment 
of the simplest model for the prey-predator interaction. They are based on 
the Adomian descomposition method [2], the time-series expansion [3], the 
homotopy-perturbation method [4, 5] , and the variational-iteration method 
[6]. We have already shown that the time-power series given by the imple- 
mentation of the homotopy-perturbation method proposed by Chowdhury 
et al [5] completely fails to yield the main features of the population dy- 
namics [7]. In this short communication we briefly analyze the results of the 
other proposals in terms of what one expects from an approach designed to 
solve problems of population dynamics [1]. 
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The chosen model is the prey-predator system [2-4, 6] 



dx{t) 



x(t)[a — by(t)], a,b>0 



dt 



dy{t) 



y{t)[c — dx(t)], c,d>0 



(1) 



dt 



where x(t) and y(t) are the populations of rabits and foxes, respectively. This 
nonlinear system exhibits a saddle point at (x s ,y s ) = (0,0) and a center at 
(x s , y s ) = (c/d, a/b) [1]. Besides, the populations obey the following curve in 
the x — y plane: 



where Xq and y$ are the initial populations at time t = 0. 

The homotopy-perturbation method proposed by Rafei et al [4] leads to 
time-series expansions of the form 



that are exactly the same as those otained earlier by Biazar et al [3]. Ob- 
viously, they are suitable about the point (a;o, yo) and will not predict the 
main features of the population dynamics (revealed at a much greater time 
span) which is what really matters in this field [1] . The variational-iteration 
method [6] yields more complicated expressions but the authors state that 
they agree with the time-power series mentioned above. 

Fig. Q] shows the populations for Case I (a = b = d = 1, c = 0.1, xo = 14, 
yo = 18) in a time span larger than that considered in the earlier studies 
already mentioned above [2-4,6]. We clearly appreciate that all those ap- 
proaches predict a wrong behaviour of the populations. Besides, the time 
span in Fig. [1] is still not large enough to cover the whole portrait of the 
system evolution. 

Case V (a — b — c — d — 1, Xq — 3, yo — 2) is much more interesting [3,4]. 
Fig. [2] shows the results in the x — y plane. The time series [3,4] predict a 
completely wrong behaviour and the resulting curve cross itself which can 
never happen as it is well known [1]. 

Summarizing: present analysis clearly shows that the Adomian descom- 
position method [2], the straightforward time-series expansion [3], the homotopy- 



ln (x c y a ) -dx-by = In (x^yfi) - dx - by 



(2) 



x{t) 

y(t) 



x + x (a - by )t + . . . 
Vo + yo(dx -c)t+ ... 



(3) 
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perturbation method [4, 5] , and the variational-iteration method [6] are com- 
pletely useless for a reasonable prediction of the evolution of even the simplest 
prey-predator systems. 
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Figure 1: Exact (solid, e) and series (dashed, s) populations for Case I 




Figure 2: Exact (solid) and approximate (dashed) populations for Case V in 
the x — y plane 
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